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ABSTRACT 


A  non-linear  theory  for  the  calculation  of  the  flow  field  of  an  oblique 
flat  plate  under  blockage  condition  is  given  using  the  techniques  of  integral 
equations.  Numerical  results  are  obtained  with  the  aid  of  a  high  speed 
digital  computer  for  the  plate  situated  mid-channel  at  values  of  the  angle 
of  attack  from  5°  to  90°  and  the  channel  width-chord  ratio  from  3  to  20. 
Also  obtained  are  results  for  the  plate  situated  at  two  different  off-center 
positions  for  a  channel  width-chord  ratio  5  and  angles  of  attack  less  than 
30° 


NOMENCLATURE 


leading  edge  of  the  plate;  origin  of  the  coordinate  system. 

trailing  edge  of  the  plate;  x  =  1. 

scale  factor. 

drag  coefficient. 

lift  coefficient. 

stagnation  point  on  the  plate. 

drag  force. 

complex  potential. 

normal  distance  from  A  to  the  upper  channel  wall. 

widths  of  the  flow  above  and  below  the  dividing  streamline  at 
upstream  infinity. 

widths  of  the  flow  above  and  below  the  cavity  at  downstream 
infinity. 

flow  at  upstream  infinity. 

flow  at  downstream  infinity  between  the  upper  cavity  and 
channel  wall. 

flow  at  downstream  infinity  between  the  lower  channel  and 
cavity  wall. 

parameters  in  the  transformation, 
free  stream  pressure, 
cavity  pressure. 

uniform  up-  and  downstream  velocities  respectively. 

channel  width  =  H  +  H  . 

1  2 

complex  velocity. 

X  +  iyi  the  physical  plane, 
angle  of  attack. 

^  '  cavitation  number, 

density  of  the  liquid. 


direction  of  the  flow. 


Introduction 


The  two-dimensional  cavity  theory  of  an  unbounded  fluid,  with  the 
recently  published  works^^^*  ^^^on  the  non-linear  solution  of  bodies  of  gen¬ 
eral  shape,  can  be  considered  a  well  established  field.  For  an  experi¬ 
mentalist  who  has  to  perform  the  tests,  however,  the  existing  theories  canr 
not  always  be  applied  directly  and  have  to  be  modified  mainly  because  the 
flow  one  creates  for  experiment  is  often  bounded  by  different  types  of 
boundaries.  These  boundaries  can  be  free  surfaces  of  constant  pressure 
if  the  equipment  in  use  is  a  jet,  rigid  walls  if  it  is  a  water  tunnel,  or  both 
if  it  is  a  free  surface  channel.  The  essential  problem  is  to  determine  the 
effect  of  the  boundaries.  In  the  past,  many  papers  have  been  published 
on  this  problem.  A  small  portion  of  these  are  listed  here  as  refer- 
ences^^^'  We  shall  not  repeat  what  they  have  done  since 

these  references  are  available  and  well  known.  This  paper  also  deals  with 
this  problem,  but  the  interest  is  focused  only  on  the  unsymmetric  flow  for 
arbitrary  angle  of  attack  and  the  boundaries  considered  are  rigid  walls.  In 
a  subsequent  report,  we  shall  treat  the  cases  of  the  free  jet  and  the  free 
surface  channel. 

It  is  a  well  known  fact  that  the  cavity  length  behind  a  body  depends 
essentially  on  cr,  an  important  parameter  known  as  the  cavitation  number 
and  defined  as 

=  (P  -  P  J/i  P  U  ^  (1) 

00  c 

Here  p^  is  the  free  stream  pressure,  U  the  free  stream  velocity,  and 
p^  is  the  cavity  pressure  presumably  nearly  constant.  In  general,  the 
length  would  increase  as  a  decreases  and  with  the  presence  of  free  sur¬ 
faces,  e.  g.  in  a  jet,  o  would  become  zero  when  the  length  approaches 
infinity.  However,  this  phenomenon  is  not  observed  in  a  water  tunnel 
with  rigid  walls.  In  this  case  a  would  "each  a  finite  positive  limit  as 
the  cavity  length  increases  indefinitely  hence  cavitation  numbers  below 
are  not  attainable.  Thus  not  all  cavity  flow  conditions  can  be  modeled  in  a 
water  tunnel.  The  phenomenon  corresponding  to  a  -  a ^  in  a  water  tunnel 
is  called  blockage  or  choking  and  the  determination  of  a ^  is  therefore  one 
of  the  central  problems  in  water  tunnel  testing  of  cavity  flows. 
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In  this  paper,  the  two-dimensional  non-linear  theory  of  a  choked 
unsymmetrical  flow  over  a  flat  plate  at  an  arbitrary  angle  of  attack  is 
worked  out  using  the  techniques  of  integral  equations.  Numerical  results 
are  obtained  on  a  high  speed  digital  computer  (IBM  7094). 

Formulation  of  the  problem  and  preliminary  calculations. 

Consider  the  idealized  cavity  flow  in  a  water  tunnel  with  rigid  walls, 
depicted  in  Fig.  1.  We  assume  a  uniform  upstream  flow  with  velocity  U, 
and  a  uniform  downstream  velocity  V  as  the  cavity,  which  has  a  station¬ 
ary  interface,  approaches  its  maximum  cross-section.  The  flat  plate,  set 
at  an  arbitrary  angle  of  attack,  can  be  located  anywhere  in  the  tunnel  and 
the  stream  after  impinging  on  the  frontal  side  of  the  plate  separates  smoothly 
at  both  its  leading  and  trailing  edges.  The  plate  is  of  length  unity  or  one  may 
say  all  dimensions  are  normalized  by  the  chord.  The  coordinate  axes  are  set 
normal  and  parallel  to  the  plate  with  the  origin  at  the  leading  edge. 

If  we  call  w,  the  velocity  vector,  with  magnitude  lw|  and  direction 
6,  w  =  |w|  e  ,  then  in  our  coordinate  system,  the  uniform  velocities  at  up- 
and  downstream  infinities  are  U  e^^  and  V  e^°’  respectively. 

The  boundary  conditions  for  the  problem  are: 

6  =  a  ,  on  the  channel  walls 

6  =  0  .  on  the  plate  from  D  to  B 
6  =  It  ,  on  the  plate  from  D  to  A 
)w|  =  V,  on  the  cavity  walls. 

The  force  coefficients  at  the  choking  condition  can  be  obtained  from 
momentum  consideration  in  terms  of  the  channel  width,  the  angle  of  attack 
and  the  critical  cavitation  number; 

P  vMh^+  h^)- P  W  =  (P^-  Pj.)W-Fj5  (2) 

where  Fj^  is  the  drag,  by  conservation  of  volume,  V{h  +  Vi  )  =  WU  and 
by  definition  cr'  =  (p  -  p  )/ipU^.  Substituting 

OO  C 

Fd  =  (P^-  Pc)  W  +  pU^W-pVUW  (3) 

The  subscript  is  dropped  since  the  cavitation  number  we  refer  to  from 
now  on  is  always  the  choked  one. 
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and  the  drag  coefficient  is 


2Fr 


PU^ 


=  2w[(l  +  j)  -  /I  +  a  , 


A^,  the  area  per  unit  span,  is  equal  to  one  for  a  plate  of  unit  length. 
The  lift  coefficient  is 


Theory 


Cl=  Cj^cota. 


(4) 


(5) 


Assuming  the  flow  to  be  irrotational,  there  exists  a  complex  poten¬ 
tial  f  =  <,p  +  14*.  The  flow  field  in  the  f -plane  is  shown  in  Fig.  2.  The  two 
channel  walls  are  represented  by  the  streamlines  4^  and  -4  while  the 
dividing  streamline  coincides  with  the  negative  i^-axis.  The  stagnation 
point  D  is  chosen  as  the  origin  and  the  cut  along  the  positive  i^-axis  re¬ 
presents  the  two  branches  of  the  streamline  split  at  D.  The  leading  and 
trailing  edges  of  the  plate  are  located  somewhere  on  the  upper  and  lower 
branches  respectively.  With  reference  to  Fig.  1,  the  values  of  and 
4/^  are  given  as 


4<  =  UH  =  Vh 
1  1  1 

4>  UH  =  Vh  .  (6' 

222 

We  proceed  to  solve  the  problem  by  introducing  a  new  variable  f2, 
defined  by 

^2  -  log  —  =  T  +  i  6 
o  w 

T  =  log  ,  e  =  -arg  w.  (7  ) 

I  w  I 

The  choice  of  using  f2  as  the  dependent  variable  for  the  problem  ‘s  sug¬ 


gested  by  the  existing  boundary  conditions. 

The  flow  fields  in  the  f  and  the  i2-planes  are  to  be  connected  through 
a  parametric  variable  ^  =  t  +  in.  Consider  the  transformation 


(8a) 
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or 

£  {^)  =  [(h^+  h^)  log  (t,  -  k^)  -  log  a  -  k^)  -  log  (t  -  k^)]  .  (8b) 

The  flow  in  the  £-plane  is  mapped  into  the  upper  half  4-plane  with  the 

boundaries  on  the  real  4  or  I -axis  as  shown  in  Fig,  3.  The  upstream  I 

behaves  like  a  source  at  I  =  k  while  the  downstreams  J  and  K  be- 

o 

have  like  sinks  at  4  =  k^  and  k^  respectively.  The  net  strength  of 
these  singularities  is  certainly  zero.  The  jump  conditions  on  the  f-plane 
further  furnish  the  relations: 


and 


^  _  -uC  1 

-  'V  (k  -  k)(k^-  k  ) 

^2  "  —  (k  ■-  kj(k^  -  k^ 


Combining  Eqs.  (9)  and  (10)  yields  the  ratio 

h  k  -  k 
1  _  2  o 

r  -  r~nr  ■ 

Z  O  1 


(9) 

(10) 


(11) 


The  constant  C  is  a  scale  factor  to  be  determined  later  in  the  theory, 
has  the  following  values  along  the  4 -axis  and  0  =  0^. 

Im  =  It ,  4  <  - 1 

Im  =  0,  4  >  1 

Ren=0,  -l<|<k,k<4<l 

1  2 

Im  i2=a,k<4<k.  (12) 

1  2 

We  further  define 


i2  =  a  +  ,  (13) 

where  is  associated  with  the  unbounded  fluid,  il  is  split  up  in  this 
way  to  make  the  solution  readily  obtainable.  The  term  ,  which  re¬ 
presents  the  effect  of  the  walls,  is  added  to  the  unbounded  fluid  trans- 
cendentally  to  maintain  the  total  exact  solution.  When  the  channel  width- 
chord  ratio  is  very  large,  then  can  be  considered  merely  as  a  per¬ 
turbation  of  the  unbounded  flow. 
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With  reference  to  Figures  4a  and  4b,  has  the  form 

C08h'‘4  =  log  a  +  fi-  1)  . 

Along  the  4-axis, 

Im  =  IT,  4  <  -1 
o 

Im  0 ,  4  >  1 

Re  0,  -1<4<  1. 

For  |4|<1,  is  given  by 

Im  2^=  P(4)  =  arctan  (/l  - 


(14) 


(15) 


(16) 


The  remaining  part, 
channel  walls,  then  takes  the  boundary  values; 

4<  -1 


2^  which  is  caused  by  the  presence  of  the 


Im  2  =0, 

1 

Im  2=0, 

1 

Re  2  =  0, 

1 

Im  2  =  a  - 
1 


4>  1 

-1<4<  k 


k  <4<  1 
2 


P(4), 


k  <4<k  . 
1  2 


The  boundary  condition  for  2^  on  the  4-axis  is  shown  in  Fig.  5.  Three 

branch  cuts  appear  on  the  4-axis;  these  are  from  4=-'JOto4  =  -l,  from 

4=1  to  4  =  +“^  and  from  4  =  k  to  4  =  k  .  The  first  two  cuts  repre- 

I  2 

sent  the  plate  and  the  last  one  the  channel  walls  in  the  z-plane. 

We  now  have  a  boundary  value  problem  involving  the  unknown 
2^  (4),  which  c 
By  analytic  continuation 


2  (4)  =  -2  (4). 

The  region  is  extended  to  the  lower  half  4"Plane  with  the  Re  2^  odd  and 

the  Im  2  even  with  respect  to  the  4-axis.  We  now  introduce  the  auxi- 
1 


liary  function 


H(4)  = 


1 


v/(;-k  )(^“k^)(4"-  1) 


(17) 


function,  2^  (4),  which  can  be  formulated  and  solved  as  a  Hilbert  problem. 


(18) 


(19) 
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which  has  the  following  properties: 

i)  H(4)  has  the  proper  branch  cuts  (^  =  -1,  k  ,  k  ,  and  +1  are  branch 

1  2 

points ). 

ii)  On  the  Re  C,  axis: 


Im  H  =  0  , 

Re  H  =  0, 

Re  H  =  -  ^  - . .  , 

y/(l-k  )(k  -l)(l  -1^) 

Re  H  =  0  , 

Im  H  =  0, 


-uo  <  4  <  - 1 
-l<|<k 

1 

Im  H  =  0,  k  <|<  k 
1  2 

k  <S<1 
2 

1  <^<  (JO 


A  new  function 

G{4)  =  njc)  m;) 

is  then  formed  which  has  on  the  Re  ^  axis; 


Im  G 
Im  G 


Im  G 


0, 

_ -  P(§) 

\/(&-k  )(k-6)(l-a  ’ 

0, 


-oo<  5  <  k 

» 

k  <k 
1  1  2 

k  <i<+-c 
2 


(7 ) 

With  the  aid  of  Plemelj's  formula'  ,  the  solution  of  the  Hilbert  problem  is 
given  immediately  by 


«  (4)  =  -  i\/(4-k^)(4-k^)(4"-l) 


^  La-P(t)]dt 
(t-C)  \/(t-k^)(k^-t)(l-t2) 


(20) 


However,  the  constants  k  and  k  cannot  be  arbitrarily  chosen  for  a  fixed 
a  because  of  a  physical  condition  that  must  be  imposed  on  the  solution. 
Specifically,  the  solution  S2  must  behave  locally  like  a  stagnation  flow 
at  the  stagnation  point.  In  the  ^-plane,  the  stagnation  point  corresponds 
to  4  =  The  local  stagnation  flow  can  be  shown  to  behave  like  log  4  as 
4  approaches  infinity;  this  behavior  is  already  incorporated  into  the 
function  We  require,  therefore  that  be  finite  at  infinity.  If  Eq. 

(20)  is  expanded  in  inverse  powers  of  this  condition  is  seen  to  be 
given  by  the  requirement  that 
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^  -  P(t) 


,  \/(t-k  )(k 

i  *  2 


dt  =  0 


(21) 


Thus  for  a  given  value  of  a,  Eq.  (21)  determines  the  relation  between 

k  and  k  . 

1  2 

We  are  now  able  to  determine  a  from  fJ  in  Eq.  (20).  At  up¬ 
stream  infinity  in  the  physical  plane,  iw|  =  U  and  4  =  k^. 

*k^  (t-k^)\/(t-k^)(k^-t)(l-t2) 

From  our  previous  discussion,  a  =  V^/u^-l,  therefore  Eq.  (22)  enables 
us  to  calculate  o  when  the  integral  on  the  right  hand  side  is  evaluated. 

In  order  to  complete  the  solution,  we  shall  relate  the  geometry 
of  the  flow  to  the  undetermined  constants.  Utilizing  the  definition  of  the 
complex  velocity,  we  have  the  relation 


dz 


w 


It 


1  «(Udf 

T  ®  ar 


With  the  aid  of  Eqs.  (8a),  (13)  and  (14)  and  upon  integration. 


c 

j  V 


(;  -k  )(;  -k  )(;  -  k  ) 
o  1  2 


When  we  apply  Eq.  (24)  along  the  plate  from  A  to  B,  the  normalized 
chord  is  obtained; 


(23) 


(24) 


'B 


Za=  1  = 


C 


X 


5  + 1 5^-1 


(6 


■k^)(S  -k  )(t. 


■k  ) 
2 


«l(^) 


di 


\/¥~- 


(6 


■k  )(6  -k  )(^  -k  )  ® 
O  1  2 


(25) 


The  first  integral  on  the  right  hand  side  represents  the  distance  from  the 
stagnation  point  to  the  trailing  edge,  we  therefore  find  the  location  of  the 
stagnation  point. 

Ha»  the  height  of  A  below  the  upper  channel  wall  can  be  deter¬ 
mined  in  a  similar  way.  With  reference  again  to  Fig.  1, 
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or 


+  Im[e"“'{zj  -  z^)] 

+  V  Imj^e  (§-k^)(&  -k  )(i  -k^)  ^ 

k 

C  f  ‘  sin(ei  +  p  -  g) 

^J.  (4  -k  HI  -k  )(l -k  ) 

“1  O  1  2 


(26) 


since  ^2^  is  purely  imaginary  and  equal  to  i6^  in  the  interval  -1^4^  k^ . 
The  channel  width  W  is  expressed  through  continuity  as 


Final  Results. 


W 


(h  +  h  ) 
1  2 


V  _  (N " 


(27) 


We  have  applied  the  theory  to  calculate  the  case  of  a  flat  plate 
positioned  mid-channel  for  various  values  of  a  and  W.  The  range  in 
a  is  from  5°  to  90°  and  in  W  from  3  chords  to  20  chords.  We  also 
obtained  results  for  two  off-center  positions  of  the  plate,  =  1.  0 
and  4.  0  for  W  =  5  chords  and  values  of  a  up  to  30°. 

The  lift  coefficient,  drag  coefficient  and  choking  cavitation  num¬ 
ber  are  plotted  versus  a  for  W  =  5  chords  and  the  plate  in  mid¬ 
channel  in  Fig.  6.  The  dash-dot  curve  in  Fig.  6  is  the  lift  coefficient 
given  by  Wu's  unbounded  fluid  theory^^\ 

In  Fig.  7,  we  show  the  enlarged  portion  of  Fig.  6  for  a  <  30° 

(8) 

with  experimental  values  from  Wade'  '  corresponding  to  the  choking 
o  determined  by  the  present  theory.  Wade's  experiment  was  performed 
with  and  the  experimental  cavitation  numbers  were  obtained  from 

measured  cavity  pressure.  Since  it  is  difficult  to  operate  the  tunnel  near 
choking  condition,  Wade  was  not  able  to  obtain  data  at  such  low  cavita¬ 
tion  numbers:  Therefore,  we  show  values  extrapolated  from  his  data 
curves.  One  notices  that  all  points  are  enveloped  by  the  two  theoretical 


curves. 
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/K) 

Cohen,  Sutherland  and  Tu'  '  have  worked  out  a  linearized  theory 
for  a  flat  plate  at  small  angles  of  attack  and  calculated  o  for  mid-channel 
position.  However,  their  definition  of  mid-channel  is  different  from  ours 
as  we  define  mid-channel  when  the  center  of  the  plate  lies  on  the  center - 
line  while  they  locate  the  leading  edge  of  the  plate  on  the  centerline.  As  a 
result,  it  is  not  possible  to  make  a  direct  comparison  between  the  two 
theories.  For  very  small  angles  of  attack,  however,  the  discrepancy 
caused  by  the  difference  in  definition  becomes  small  and  for  a  =  6°,  W  =  5, 
the  present  theory  gives  O-  0.133,  0.198  for  0.959 

and  their  theory  gives  o  ~  0.  1 3  and  22  for  H^/(W-H^)=  1.  The 

values  of  their  calculations  are  read  from  Fig.  10  of  Ref.  3,  hence  are 
only  approximate. 

Figure  8  shows  the  effect  of  plate  position  for  W  =  5  chords.  In 
this  figure  C^,  and  o  are  plotted  versus  a  for  a  <  30°  and  two 

values  of  H^.  The  dashed  curves  are  =1.0  and  the  solid  ones  are 

Ha  =  4.  0.  As  would  be  expected  Cj^,  and  c  all  increase  as  the 

plate  is  lowered  in  the  channel. 

We  show  C.  and  o  versus  a  for  W  =  8  with  experimental  values 

/Q\  ^ 

from  Parkin'  '  in  Fig.  9.  Again  the  experimental  cavitation  numbers  were 
obtained  from  measured  cavity  pressure,  however,  in  this  case  only  the 
data  for  a  =  8°  and  10°  were  obtained  from  extrapolated  curves. 

Figures  10a  and  b  show  drag  coefficient,  and  choking  cavita¬ 

tion  number,  cr,  versus  channel  width,  W,  for  a  flat  plate  located  mid¬ 
channel.  The  range  in  angle  of  attack  is  5°  to  90°  and  in  W  from  3  to  20 
chords.  C.  may  easily  be  obtained  from  Fig.  10a  and  Eq.  (5)  and  there- 
fore  is  not  presented  here. 

In  order  to  show  the  effect  of  the  channel  walls,  we  have  plotted 
RCj^,  the  ratio  of  bounded  to  unbounded  (RCj^  =  RCj^  = 

versus  the  channel  width  W  in  Fig.  11.  The  curves  are  for  the  same 
values  of  a  and  W  as  Figs.  10a  and  b  except  that  for  a  =  90°  RCj^  is 
shown  for  W  as  low  as  1.  224.  As  would  be  expected  when  W  =  20  there 
is  very  little  wall  effect  and  as  the  channel  becomes  smaller  the  ratio 
drops  until  the  walls  approach  the  plate  where  the  ratio  increases.  This 
effect  is  shown  for  P  =  90°  where  the  ratio  tends  to  infinity  at  W  =  1. 
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The  curve  of  RCj^  for  each  a  will  have  a  similar  asymptote  at  W  =  sin  a. 

In  Fig.  12,  we  have  shown  the  distance  of  the  stagnation  point  from 
the  leading  edge  as  a  function  of  a.  At  a  =  15°,  it  is  already  difficult  to 
tell  the  difference  between  the  two  points  and  at  a  =  ^  the  stagnation  point 
is  at  mid-plate  as  one  would  expect. 
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APPENDIX 

Procedures  for  numerical  computations 

In  a  direct  problem  one  would  specify  the  channel  width,  W,  and  the 

location  of  the  plate  (of  chord  unity),  H^,  for  a  given  angle  of  attack,  then 

calculate  o,  h  ,  h  and  the  force  coefficients.  However,  in  view  of  the 
1  2 

form  of  our  solution  it  is  necessary  to  do  the  inverse  problem  of  choosing 

the  transformation  parameters  k  ,  ,  and  k^  (the  downstream  conditions, 

h  and  h  ,  when  the  flow  is  choked)  and  calculating  W  and  H .  .  The 
12 

values  of  h  and  h  are  related  to  the  k's  through  the  transformations 
in  Eqs.  (8a),  (8b)  and  the  jump  conditions  in  Eqs.  (9)  and  (10).  The 
choice  of  k  defines  the  geometry  so  that  the  scaling  factor,  C,  or  more 
precisely,  the  ratio  V/C  is  determined  from  Eq.  (25). 

Given  a  value  of  k  ,  k^  is  easily  determined  by  iteration  of  the 
integral  condition  of  Eq.  (21).  The  integral  in  this  equation  cannot  be  num¬ 
erically  integrated  in  its  present  form  so  we  divide  it  into  two  integrals 


p.  _ 


a-i3(t) 


-k  +k 

r_i _ 2 

dt  =  \  2 


('k 

dt.  \  ^ 


a-P(t) 


,  „■  V-.- - ut  =  \  2  T  .  QX  r  ^  —  '7— = 

)(k-t)(l-t^)  ^(t-k^  )(k^-t)(l-tM  ‘^^+k^  i/(t-k^)(k^-t)(l-t2) 

2 


dt: 


make  the  transformations 

p  =  +  k. 


s  =  k  - 


in  the  first  integral, 
in  the  second  integral. 


and 


x„"  =  i(k  -  k  ) 


2  1 


and  Eq.  (21)  becomes 


p-o  d, ♦  f°  d,  =  0. 

‘^o  p^^Zx^-  X*  •'o  s^^Zx^-  x^ 


(A-1) 


(A-2) 


This  integral  iteration  presents  no  problem  since  all  calculations  are  pre 
formed  on  a  high  speed  digital  computer  (IBM  7094).  We  now  can  calcu¬ 
late  the  remaining  transformation  parameter,  k^,  from  the  specified 

ratio  h  /h  and  Eq.  (11). 

2  1 
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To  aid  in  the  trial  and  error  solution  of  this  problem  the  behavior 
of  and  for  the  plate  in  mid-channel  are  shown  in  Fig.  A-1  and 
for  the  off  center  positions  in  Fig.  A-2.  The  remaining  parameter  to  be 
guessed  is  h  /h  which  may  be  approximated  by  (W-H. -sina  )/H. .  For 
the  mid-channel  case  this  quantity  is  roughly  unity. 

The  choking  cavitation  number,  a,  the  channel  width,  W,  the  drag 
coefficient,  Cj^,  and  the  lift  coefficient,  are  readily  determined 

from  Eqs.  (22),  (27),  (4)  and  (5)  respectively.  The  integral  in  Eq.  (22) 
is  treated  in  the  same  way  as  Eq.  (21). 

To  determine  we  must  first  evaluate  the  ratio  V/C.  From 

Eq.  (25) 


V 


00 

=  "It 


)(4 


k  )ii  -k  ) 
1  2 


■I 


00 


4+l/P^l 

(t+k  )(Uk  )(4+k  ) 

o  1  2 


n'd) 

e  1  d5 , 


(A-3) 


where  (4)  is  given  in  Eq.  (20)  and 


«/(!)=  -  7)/(5+k  )(4  +  k  )(|2.  1) 


g  -P(t) 


(t+e)V(t-k  )(k  -  t)(i  -t") 

1  2 


dt. 


(A-4) 


For  numerical  integration,  it  is  more  convenient  to  introduce  the  new  in¬ 
tegration  variable  t  =  4  The  first  integral  on  the  right  hand  side  of 
Eq.  (A-3)  now  may  be  written 


00 


(4  -k  )(4  -k  )(4  -k  ) 

O  1  2 


e“<«> 


(*  °  1  +  \/l  -T^  ^  ^  f  1  +  \ll  -T^  i?  “  J 

J  n-J^HT^w-iFr)  ^  ‘  '^^jri-FT)(r-kT)(i-kT)«^ 

O  °  1  2  J  O  1  2 


(A-5) 


and  the  second  integral 
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r  °  1 + \/i  -  n .  r  1  +  \/i  - 

i(i  +  kTj(i  +  kT)(i  +  kTj  ®  +  J  (irk;T){i+k-TKr+kT)- 

OOi2  ^Oi2 


dr. 


(A-6) 


where 


a 


=  ^\/(l  -  k^T)(l  -  k  T)(l  -  T^)  [g^+  Tg^+  T2g^+  ...+  T^g^j 


(A-7) 


1  \I{1  +  k  t)(1  +  k  t)(1  -t*) 

w  1  2 


-  Tg  +  T^g 
2  3 


+  T 


(A-8) 


_ a  "  P(t) _ 

(1  -tT)\/(t-k  )(k  -  t)(l  -t^) 
12 


dt, 


{A-9) 


~\/(l  +  kT){l+kT)(l-r‘)  f 


-P(t) 


4  (l  +  tT)\/(t-k  ){k  -  t){l  -t*) 


dt, 


and 


=  t"[a-P(t)Jdt 
"  {y^t-k)(k-t)(l  -  t^) 

^  I  Z 


(A-IO] 


(A-11) 


and  iJ  ^  are  obtained  from  the  power  series  expansion  of  T  and  there¬ 
fore  are  suitable  for  the  calculation  of  small  values  of  t  or  large  values  of 
For  a  properly  chosen  T^,  say  t^=  0.  1,  a  truncated  series  of  seven 
terms  is  sufficient  to  give  (t^)  identical  to  to  five  significant 

figures. 

The  remaining  quantity  to  be  evaluated  is  the  integral  in  Eq.  (26) 
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i 

I 


sin  (6  +  p  -  a ) 

.  (s-kJ(Uk  )ft-k  ) 


Some  analysis  is  needed  before  a  straightforward  numerical  integration  can 
be  performed.  One  notices  that  the  denominator  of  the  integrand  has  a 
factor  (I  -  k^)  which  causes  the  integral  to  be  logarithmically  singular  un¬ 
less  the  numerator  also  vanishes  at  the  upper  limit.  It  is  therefore  neces¬ 
sary  to  examine  the  behavior  of  the  numerator  as  4  approaches  k  .  Let  us 
study  the  behavior  of 


_ (.1 

e  =  ^  7(k  -  4)(k  -  4)(l  -4")  j 


Q  -  P(t) 


(t-4)\/{t-k  Kk^-  1^<1  -t^) 


dt 


(A-12) 


as  4  approaches  k  .  We  denote  the  regular  part  of  the  integrand  in  the 


neighborhood  of  t  =  k  as 

1 


F(t)  = 


a-P(t) 


\/{k^-t)(l  -  t^) 


(A-13) 


F(t)  has  the  Taylor  expansion 


,  F''(k  )  2 

F(t)  ^  F(k  )  +  F  (k  )(t-k  )  +  --Mt-k)  + 

i  11^  i 


where  the  coefficients  F(k  ),  F*(k  )  etc.  are  finite  numbers  which  can  be 

1  1 

evaluated  without  any  difficulty.  If  we  now  subtract  from  F{t)  the  first 

terms  of  its  Taylor  expansion  then  add  them  on  later,  the  integral  in  6 

can  be  written 
k 

"  [a-P(t)Jdt _ 


I 


(t-4)\/(t-k  )(k^-t)(l-t2) 


dt 


(t-4)\/Fir 


F’'(k  ) 

F(t)-F(k  )-F'(k  )(t-k  ) - 2-i- 

1  11^ 


(t-k  ) 

i 


k  I  k 

P  (t-k  )"2  r 

-TT-'**  — 


(t-k  ) 


‘  dt  +  - 2  ^ 


k  j 

F"(k  )  C*  (t-k  )2 


I 


IT 


dt. 


(A- 14 
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The  integrand  of  the  first  integral  of  Eq.  (A-14)  is  finite  and  tends  to  zero 
as  t  approaches  k  because  the  terms  in  the  bracket  are  of  the  order 

3  1 

(t-k  )  .  The  remaining  integrals  of  Eq.  (A-14)  are  integrated  in  closed 
1 

form: 


r  JL  dt 

J  t  -I  \ft  -  k 


k  -§ 

r  -  \ 

i 

k  -  k 

L  V 

2  1 

i 

k 


\^t-k 

j  t^r  dt  =  2i/k^-k  -  2|/k  -  I 


X 

1 

_ r 

r  -  tan  S 

i 

k  -k 

\ 

2  1'  _ 

k  3_ 

'  ^  (t-k  )2 


-4-  dt  =  i  (k  -  k  )2-  2(k  -k  )l/k  -  ^  +  2(k 

-b  J2l  211  1 


” 

k  -  e 

IT  -1 

j-  -un  , 

i 

k  -  k 

1 

2  1  _ 

From  Eqs.  (A-13)  and  (A-14)  it  can  be  shown  that  in  the  neighborhood  of 
§  =  k  ,  6  has  the  asymptotic  expression 

6  =  a  -  (3  +  MJ)fi  -  i  +  0(k  -  &) 

1  ’1  1 


where 


M  = 


V/(k  -  k  )(1  -k/) 


dt 


TT 


k 

1 


F(t)  -  F(k^) 


F'(k  )(t 


"(k)  2- 

^(t-k) 


2  ^ 

^  Jk  -  k 

^  2  1 

In  our  final  form  for  numerical  computation,  the  range  of  integra¬ 
tion  of  the  integral  in  Eq.  (26)  is  split  into  two  intervals; 


2F'(k) 


\/l  -  k 


(k  -  k  )  + 
1  2  1 


F"(k)  ^ _  2 

^5— i-  Jl  -  k  ^  (k  -  k  ) 
3tt  w  1  '  2  1 


(A-15) 
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k-« 

sin  (6  +  P  -  a )  Z'  *  si 


(i-kj(i-k)(e-k) 


sin (6  +  p  -  a) 

k  ){§-k  )(§  -k) 
O  1  2 


d|  + 


k 

S' 

k-€ 
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(6  +  p  -  a)  di 

X 

(§-k  )(§-k  ){l-k) 

O  i  2 


(A-16) 


where  «  ,  small  and  positive,  is  chosen  so  that  at  |  =  k  ,  sin(6  +  p  -  a) 

1  1 

can  be  safely  approximated  by  0^  +  P-a.  If  we  replace  6  ,  in  the  second 
integral  of  Eq.  (A-16)  by  its  asymptotic  expression,  Eq.  (26)  becomes 


H.  =  h  + 
A  1 


C 

V 


sin(6^  +  P  -  a  ) 
-k)(^-k)(l  -k) 

O  1  2 


di  -  2M 


I 


dX 

(k^-§)(k-l)  . 


(A-17) 


where 

=  k  -  ^ 

1 

The  integrals  in  Eq.  (A-17)  can  be  easily  evaluated  since  their  integrands 
are  regular  everywhere. 


Fig.  5  Boundary  values  of  SI  on  Re  4~axis. 


Fig.  7  Enlarged  section  of  Fig  6  for  small  angles  of  attack  with 
experimental  data  extrapolated  from  Wade's  data  curves, 


Fig.  8  The  effect  of  angle  of  attack  on  C,  ,  C„  and  a  for  a  flat  plate  Fig.  9  Comparison  of  C.  from  the  present  theory  with  the  unbounded 

L.  D  ^  i-' 

located  in  two  off-center  posititins  in  a  channel  of  W  =  5.  theory  and  experimental  data  from  Parkin  for  a  plate  located 

mid-channel  in  a  channel  of  W  -  8. 


Effect  of  channel  width  on  (a)  the  drag  coefficient  and  (b)  the  choking  cavitation  number  for  a  flat  plate 
located  mid'C'naiuxel  at  various  angles  of  attack. 


Fig.  11  The  ratio  of  (bounded)  to  Cp  (unbounded)  versus  W  for  a  flat  plate  located  mid-channel. 


Fig. 


Distance  of  stagnation  point  from  the  leading  edge  for  W  =  5  and  the  plate  in  mid-channel. 


PLATE  IN  MID-CHANNEL 


flat  plate  located  mid-channel. 


Fig.  A-2  Effect  of  angle  of  attack  and  plate  position  on  the  transfornnation  parameters  (a)  k  and  (b)  k  for 
a  flat  plate  in  two  off-center  positions  H.  =  1  and  4,  in  a  channel  of  W  =  5.  ^  * 
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